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Abstract 

We investigate a reaction of boron trichloride (BCI3) with iron(III) hydroxide 
(Fe(OH)3) by ab initio quantum chemical calculation as a simple model for a reaction 
of iron impurities in BCI3 gas. We also examine a reaction with water. We find 
that compounds such as Fe(Cl)(OBCl 2 )2(OHBCl 2 ) and Fe(Cl) 2 (OBCl 2 )(OHBCl 2 ) 
are formed while producing HC1 and reaction paths to them are revealed. We also 
analyze the stabilization mechanism of these paths using newly-developed interaction 
energy density derived from electronic stress tensor in the framework of the Regional 
DFT (Density Functional Theory) and Rigged QED (Quantum ElectroDynamics) . 



1 Introduction 



In this paper, we investigate a reaction of boron trichloride (BCI3) with iron(III) hydroxide 
(Fe(OH) 3 ) by ab initio quantum chemical calculation. The purpose of the paper is two fold. 
One is to present how Fe(OH) 3 reacts in BCI3 gas which could be relevant to an industrial 
process. Secondly, by applying newly-developed interaction energy density concept to the 
reaction, we would like to test its validity in particular regarding its ability to describe 
stabilization through chemical reaction. 

Let us start from describing some industrial background. Boron trichloride is one of 
the semiconductor gases, which is used, for example, as a preferential plasma etching gas 
for aluminum and a source of boron for p-type doping in the process of chemical vapor 
deposition [H El [3]. High purity is required for semiconductor gases to be used in the 
production process of recent high integrated device and thin insulating film. Required 
impurity level in recent years has been lower and lower. In the future, it is expected 
to be parts-per-trillion (ppt) level [4J. To achieve lower impurity level, it is necessary 
to remove all causes which might contaminate the gases in the whole process including 
gas transportation. As examples of concerned impurities, there are water and metal. 
Once water or metallic impurities are mixed in BCI3, they can deteriorate the product 
performance and the process yield. Generally, to remove impurities from ultra high purity 
gases, ceramic or synthetic resin filters are used. However, metallic impurities in boron 
containing gases, including BCI3 gas, are not removed well through these types of filters [5] 
and reason for this is not known. It may mean that metallic impurities in boron containing 
gases have structures which easily pass through these filters as the result of the interaction 
between gases and metal. So far, methods for removing impurities in high purity gases 
have been discussed in the literatures, but structures or states of impurities have not been 
studied. In particular, it is difficult to directly observe the structure of impurities in gas 
phase. Therefore, we consider that it is worthwhile to investigate the molecular states of 
impurities by a computational method as a precursory study toward establishing a more 
effective method for removing impurities. 

Among several metallic impurities, since most abundant one is iron, we focus on iron- 
including impurities. We can think of several possible sources for iron impurities. The 
prime suspect is rust from the welded spots in the ductwork. Although the ductwork for 
semiconductor gases has corrosion resistance, the welded spots are relatively weak against 
halogen-contained gases, especially, under the presence of water [6j. Note that at most a 
few parts-per-million (ppm) of water is mixed in the ductwork. 

As a first step to guess how iron impurities react with BCI3 gas, we consider a reaction 
of Fe(OH) 3 with BCI3. The first reason why we pick up this iron compound is that we 
would like to see a reaction with a hydroxyl function since hydrochloric (HC1) gas is known 
to exist in BCI3 gas. Such HC1 gas may come from the reaction with the iron impurities. 
The second reason is that rust mentioned above is likely to be in the form of goethite 
(a-FeO(OH)) [3 El IS] ? and its monohydrated form can be described as Fe(OH) 3 . 

This paper is organized as follows. In the next section, we briefly explain our quan- 
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turn chemical computation method. We also describe our analysis method based on the 
Regional DFT (Density Functional Theory) and the Rigged QED (Quantum ElectroDy- 
namics), and in particular we define the interaction energy density. In Sec. [3| we show 
our results on the reaction of BCI3 with water and that with Fe(OH) 3 . The final section 
is devoted to our conclusion. 

2 Calculation Methods 

2.1 Ab initio electronic structure calculation 

We perform ab initio quantum chemical calculation for several chemical reactions us- 
ing density functional theory. In this calculation, we adopt the Lee- Yang-Parr (LYP) 
[TO] gradient-corrected functionals for the correlation interaction, and it is employed with 
Becke's hybrid three parameters [11] for generalized-gradient-approximation (GGA) exchange- 
correlation functions (B3LYP). We employ 6-311G* basis set for Fe, B, O and CI, which 
consists of all electron basis set by Wachters-Hay [121 [13], an d 6-311G** basis set for H. 
This calculation shows reasonable results for high-spin states of iron hydroxide clusters. 
Analytical vibrational frequencies are obtained to calculate the zero-point energy (ZPE) 
correction. In this work, the geometric optimized structures and the electronic structures 
of each cluster model are calculated by GAUSSIAN03 program package [141. The part of 
visualization in this paper is done using MOLDEN[15] and VMD[T6] softwares. 

2.2 Interaction energy density analysis 

In the following section, we use newly-developed interaction energy density in our labora- 
tory to analyze how and in which part of molecules are (de-)stabilized during the chemical 
process. This quantity is used in Ref. [17] to describe the stabilization of molecules through 
covalent bonds and van der Waals bonds. (It is also used in Ref. [18] recently.) It is defined 
in the framework of the Regional DFT and the Rigged QED [HI EDI EH [221 ESI E! ESI EH] 
and can be calculated from the electronic stress tensor density *r* s (r) whose components 
are given by 



where {fc, 1} = {1, 2, 3}, m is the electron mass, and Vi and ipi(f) is the occupation number 
and natural orbital of the zth state, respectively. Taking a trace of *r* s (f) gives energy 
density of the quantum system at each point in space. The energy density £^(f) is given 




d 2 ipi (f ) di/>* (f) dipi (r ) <9 2 ^* (f) 
dx k dx l dx k dx l dx k dx l 
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Now, we can define the interaction energy density. Suppose that the system under 
consideration formally consists of two parts A and B and that the whole system has the 
energy distribution as s^ AB (r). When the parts A and B are considered separately, they 
have the energy distribution ef A {f) and ef B {f) respectively and sf A (r)+ef B {f) ^ ab(^0- 
The difference stems from stabilization or destabilization due to the reaction between A 
and B at each point in space and we call it the interaction energy density Asf(r). Namely, 

Aef (f) = <AB(fO - (< A (r) + < B (rl} . (3) 

The region with negative Ae^(r) corresponds to the stabilized region and the positive 
region denotes the destabilized region. 

We also use conventional electron density difference An(r) for the later discussion as 
defined below. 

An AB (f) = n AB (f) - {n A (f) + n B (r)} , (4) 

where n(r) is the ordinary electron density at r. 

As shown in Ref. [T7j, in covalent bonding, the larger stabilization accompanies the 
larger atomic population (the region with negative Asf(r) has positive An(r)). However, 
note that this is not the case for van der Waals bonds, in which the stabilized region 
corresponds to the decrease in the electron density [17J (the region with negative Asf(f) 
has negative An(r)). 

We use Molecular Regional DFT (MRDFT) package [27J to compute these quantities. 

3 Results and discussion 
3.1 Reaction with water 

In this section, we study the reaction of BCI3 with water. Although the hydrolysis of BCI3, 
BCl3(g) + 3H 2 0(1) — ► B(OH) 3 (aq) + 3HCl(aq), is a well-known textbook-level reaction, 
the situation we consider here is slightly different. In our case, since H 2 is much less 
(ppm level) than BC1 3 , a relevant reaction would be BC1 3 + H 2 BCl 2 (OH) + HC1. 

Then, we start by examining how H 2 approach BCI3 and form a complex CI3B — OH 2 . 
We take the distance between B and O as a parameter, D, and for several values of D, we 
calculate the optimized configuration of the other atoms. 

The obtained structures are shown in Fig. [T] We can regard them to be the snapshots 
of the continuous reaction process. While D is between 4.8 A and 4.0 A, we find that H 2 
approach with a small angle from the BCI3 plane. In Fig. [2j we show relative energy and 
charge transfer as functions of D. The charge transfer is calculated from the Mulliken 
charge. We obtain the structure of the stable reactant complex CI3B — OH 2 as shown in 
Fig. [3] at D = 1.693 A. Below, we refer to this structure as "RC". This geometry is in 
good agreement with the one in the literature [28l [29] . We find that when H 2 approaches 
BCI3 from an infinite distance, there is no energy barrier and total energy is stabilized 
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by 0.2607 eV. Also note that the reaction proceeds as charge is transferred from BCI3 to 
water (electrons from water to BCI3). 

Next, we search a reaction path from RC until the detachment of HC1. The energy 
along a certain intrinsic reaction coordinate (IRC) is plotted in Fig. [3J Some intermediate 
structures are also shown. In particular, the one labeled 1 is RC and 3 is the transition 
state (TS) respectively. This is reorganized in Fig. [5] and Table [l] which show the reaction 
pathway and relative energy. The activation energy is found to be 0.7300 eV. By detaching 
HC1, the system stabilizes by 0.6799 eV from RC. In the final step, there is energy increment 
of 0.0990 eV which corresponds to the strength of the hydrogen bond between H in HC1 
and O in BCl2(OH) (the distance between H and O is 2.054 A). However, since it is rather 
small, BCl2(OH) and HC1 are expected to be separated away in the gas phase. 

Now, we analyze this reaction by using interaction energy density Aef(f) and electron 
density difference An(r) as introduced in Sec. [2| We first examine the process from 
BCI3 + H 2 to RC, namely the reactant complex formation. This is shown in Fig. [6] From 
An{f) of panels (a) and (b), we see that BCI3 and H 2 are polarized by the existence of 
the other molecule when they are separated. Then they attract each other by electrostatic 
interaction. After they approach closely as panel (c), electrons move from H 2 to BCI3 
rapidly as shown by the expansion of pink regions around BCI3 in panels (d), (e) and (f). 
This is consistent with the charge transfer plot in Fig. [2| which shows steepening trend for 
D < 3.5 A. As for Aef (r), general feature we notice is that positive An(r) region (colored 
in pink) roughly corresponds to negative Asf(r) region (colored in blue) and vice versa. 
This indicates that each part in the system is stabilized by the increase in the electron 
density just as in the usual covalent bonding. 

The process from RC to BCl 2 (OH) • HC1 via TS is next analyzed and results are shown 
m Fig. We again see the correspondence between positive An{r) region and negative 
As^(r) region. Since we partition the system into BCl2(OH) and HC1 for calculating 
An(f) and Aef(r), it is easier to see the process backward from the panel (g). We see 
that red destabilized region expands from panel (g) to (a). This destabilized region is 
especially large in TS (panel (c)) around H 2 0. Although this partitioning is not so well 
defined in (a) and (b), we see destabilized region around H 2 too. This is consistent with 
the energy level relation that RC has higher energy than BCl 2 (OH) • HC1. 

3.2 Reaction with iron hydroxide 

In this section, we study the reaction of BCI3 with Fe(OH) 3 . Fe(OH) 3 is our model of 
iron impurity and we investigate how this can react with BCI3 molecules to produce HC1. 
It should be mentioned that Fe(OH) 3 has lower energy than FeO(OH)-H 2 0. As shown in 
Table [2| since sextet is the most stable, we adopt this state in the following calculation. 

The entire reaction path we have searched is shown in Fig. [8) Table [3] shows relative 
energy for each step. For convenience, we split this path into four as in Fig. [9] (a)-(d) 
showing structures for each step. We will give detailed description for each of them below. 

The earliest stages in the reaction path are shown in Fig. [9] (a). First, it is shown that 
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BCI3 approaches Fe(OH) 3 with some energy stabilization forming a reaction complex S2. 
A possible simple path to detach HC1 from S2 is shown as S2 —> TSi — > S3. This is just 



like the one we have seen in the reaction between BCI3 and H 2 described in Sec. |3.1[ We 
need the activation energy of 1.1653 eV. In the case of BCI3 and Fe(OH) 3 , however, the 
reaction complex S2 can turn into more stable structure by cutting a bond B — CI in the 
ring consists of B, CI, O, and Fe atoms and leaving a Fe — CI bond (S2 — > S4). Since the 
activation energy is rather high, it is likely that the latter path is taken in the BCI3 gas. 
Also, as we will see below, the presence of many BCI3 molecules (relative to Fe(OH) 3 ) 
makes it possible to open more effective paths to detach HC1 molecules. 

Now, let us consider the case that another BCI3 approaches S4. One more bond between 
BCI3 and OH is formed to make S5 and energy is stabilized by 1.1178 eV. To detach HC1 
from S5, there are two reaction paths as shown in Fig. [9] (b). One is from S5 to S6 
via TS2a and another is via TS2b- The former is similar to the BCI3+H2O case or the 
path S2 —> TSi — > S3 mentioned just above and its activation energy is relatively high, 
0.8665 eV. By contrast, the latter has much lower barrier of 0.15083 eV. In this path, since 
HC1 is detached from CI and OH which are connected to different B, the distance between 
H and CI is shorter in TS and much less energy is needed to form the bond. 

The path after that is opened in a similar way (Fig. [9] (c)). Namely, HC1 is not likely to 
be produced from within S7 but detach after a complex with one more BCI3 is formed. This 
complex is shown as S$. Then HC1 is again formed from CI and OH which are connected 
to different B with a relatively low activation energy of 0.08919 eV (Ss — > TS3 — ► S9). 

So far, we have learned that we can find an energetically more favorable path to detach 
HC1 by attaching BCI3 beforehand. Then, similar path to S2 —> S4 — > S5 may take place 
for S5 by attaching BCI3. This turns out to be true and we find a path S5 —> Sn — > S12 
as shown in Fig. [9] (d). Sn is stabilized by 0.59024 eV from S5 by forming another Fe — CI 
bond. After that, S12 is stabilized by 0.8491 eV from Sn by the interaction between BCI3 
and the last OH group bonded to Fe. 

Finally, how HC1 is detached from S12 is described. As shown by the path S12 — ^ TS4 — > S13, 
CI and OH which are connected to different B bond to form HC1. This is similar to 
S5 — ► TS2b — ► S6 and Ss —> TS3 — > S9. It also has a relatively low activation energy of 
0.22614 eV. We further examine whether one more HC1 can be detached from S13. How- 
ever, we cannot find such a path and we find that BCI2OH is detached instead. Thus, the 
final product is Fe(Cl) 2 (OBCl 2 )(OHBCl 2 ) shown as S M . 

To sum up, we have found 

Fe(OH) 3 + BCI3 Fe(OH) 2 (OBCl 2 ) + HC1 (5) 

has high activation energy and unlikely to occur but there are two reactions with lower en- 
ergy barriers which produce HC1 and more stable Fe-compounds. They can be summarized 

as 

Fe(OH) 3 + 3BC1 3 -> FeCl(OBCl 2 )2(OHBCl 2 ) + 2HC1 (6) 
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and 



Fe(OH) 3 + 3BC1 3 -> FeCl 2 (OBCl 2 )(OHBCl 2 ) + HC1 + BC1 2 0H 



(7) 



In passing, it may be useful to comment on the geometrical structure of each complex 
in the reaction path. Generally speaking, a four-coordinate complex forms a square-planar 
or tetrahedral structure. We calculate the skewness of the complex defined as follows [30] : 



a 



opt 



Vr 



(8) 



opt 



where V is the volume of the tetrahedron defined by the four atoms directly connected to 
Fe, and V opt is the volume of a regular tetrahedron which has the common circumsphere 
to that tetrahedron. If a = 0, the tetrahedron is regular and if a = 1, it is square planar. 
The result is shown in Table [3| This result shows that the complexes we have dealt with 
are very close to regular tetrahedrons. 

The analyses using An(r) and Ae^(r) for some parts of the reaction path are shown 



in Figs. 10 12, As is the cases which are examined in Sec. |3.1[ positive An(r) region 
corresponds to negative Aef(r) region in general. 



Fig. 



10 shows the process of BCI3 
approaching a complex Fe(Cl)(OH) 2 (BCl20H) to form a reaction complex and Fig. 11 
shows the process of detaching HC1 from the reaction complex. They are respectively 
similar to what we have seen in Figs. [6] and [7j Namely, we see blue stabilized region grows 
between BCI3 and Fe(Cl)(OH) 2 (BCl20H) as the process proceeds from (a) to (e) in Fig. 10 
and red destabilized region develops in the direction from (e) to (a) in Fig. 11, We show 



in Fig. 12 one more example of As^(r) for the process of detaching HC1 from the reaction 
complex. When we look Fig. [12] in the direction from panel (f) to (a), we see that red 



destabilized region grows, as is the case of Fig. 11 



4 Conclusion 

We have investigated a reaction of BCI3 with Fe(OH) 3 by ab initio quantum chemical cal- 
culation as one of the simplest models for a reaction of iron impurities in BCI3 gas. We have 
found that compounds such as Fe(Cl)(OBCl 2 ) 2 (OHBCl 2 ) and Fe(Cl) 2 (OBCl 2 )(OHBCl 2 ) 
are formed while producing HC1. The reaction paths to them are examined in detail and 
their activation energy is found to be relatively low due to the formation of a Fe-complex 
coordinated by several BCI3 before detaching HC1. We have also examined a reaction with 
a single H 2 molecule (remember that H 2 is rare in the BCI3 gas) and have found that 
it has high activation energy. Such difference in energy barriers indicates that it is more 
likely that the observed HC1 originates from the reaction of BCI3 with iron impurities 
rather than from the reaction with H 2 0. 

We have also analyzed the stabilization mechanism of these paths using newly-developed 
interaction energy density Asf (f) in our laboratory derived from electronic stress tensor 
in the framework of the Regional DFT and Rigged QED. We have compared this with 
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electron density difference An(f). We have found correspondence between positive (neg- 
ative) An(f) region and negative (positive) Ae^(r) region in generai. This indicates a 
covaient bond that a bond is stabilized by the increase in electron density. We believe this 
interaction energy density is very usefuf to analyze and visualize how and in which part 
of molecules are (de-) stabilized during the chemical process. Integrating the interaction 
energy density over some region would give good quantitative measure of stabilization. 
This will be investigated in our future work. 

Although it is too early to conclude that the reaction paths we have shown are realized 
in the BCI3 gas in the ductwork, it is reasonable to imagine iron impurities play some role 
in producing HC1. More detailed modeling of iron impurities in future would give us more 
hints for this issue. 
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Table 1: Relative energy along the pathway of the reaction of BCI3 with water shown in 
Fig.§ 





Relative energy (eV) 


BC1 3 +H 2 





RC 


-0.26065 


TS 


0.36501 


BCl 2 (OH) • HC1 


-0.94056 


BCl 2 (OH) + HC1 


-0.84152 



Table 2: Relative energy in units of eV for the optimized structure of several spin states 
of Fe(OH) 3 and FeO(OH)-H 2 0. We take sextet of Fe(OH) 3 which has the lowest energy 
as the reference point. 





Quartet 


Sextet 


Octet 


Fe(OH) 3 
FeO(OH)-H 2 


0.01366 
0.04452 




0.04936 


5.48862 
0.17616 
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Table 3: Relative energy and skewness of tetrahedron (Eq. ([8])) along the pathway of the 
reaction of BCI3 with Fe(OH) 3 shown in Fig. [8j 





Relative energy (eV) 


Skewness of tetrahedron o 


Si 





(3-coordinate) 


s 2 


-1.83438 


0.0990 


TSi 


-0.66899 


0.0990 


S 3 


-1.70663 


0.1013 


S 4 


-2.36146 


0.0961 


s 5 


-3.47929 


0.0324 


TS 2a 


-2.61284 


0.0570 


TS2b 


-3.32846 


0.0355 


S6 


-3.54107 


0.0291 


s 7 


-3.41744 


0.0801 


s 8 


-4.27892 


0.1591 


TS 3 


-4.18973 


0.0723 


S9 


-4.57131 


0.0515 


S10 


-4.47948 


0.0549 


Sn 


-3.88924 


(5-coordinate) 


S12 


-4.73834 


(5-coordinate) 


TS 4 


-4.51220 


(5-coordinate) 


S13 


-4.72870 


0.0364 


S14 


-4.49319 


0.0492 
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t 




D=2.4A^ D=2.2A 




D=2.0A 



Figure 1: Optimized configurations of BCI3 and H 2 at the given B — O distance (D). 
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-0.4 1 1 1 1 1 1 1 1 -0.5 

1.5 2 2.5 3 3.5 4 4.5 5 

D(A) 

Figure 2: The plots of relative energy (solid red line) and charge transfer (dashed green 
line) vs B — O distance. The reference energy is the energy of structure in which H 2 
locates at an infinite distance from BCI3. The change in the charge of BCI3 is plotted as 
the charge transfer on the right axis. Namely, the negative charge transfer indicates that 
electrons moves in to BCI3 from H 2 0. 




Figure 3: The structure of the stable reactant complex (RC), CI3B — OH 2 . 
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- o L 5 \ 









Reaction coordinate 



Figure 4: Relative energy of the path from C1 3 B— OH 2 (RC) to BC1 2 (0H) • HC1. The 
reference energy is the energy of structure in which H 2 locates at an infinite distance 
from BCI3. 





TS 


BC1 3 + H 2 










\ BCl 2 (OH) + HC1 




BCl 2 (OH) • HC1 


Reaction steps 



Figure 5: Reaction pathway from BC1 3 + H 2 to BCl 2 (OH) + HC1. 
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Figure 6: Optimized configurations (left column), electron density difference An(f) (mid- 
dle column) and interaction energy density Aef(f) (right column) during the process of 
the CI3B — OH 2 (RC) formation. We partition the system into BCI3 and H 2 for calculat- 
ing An(r) and Aef(r). For An(r), the light blue region has negative value and the pink 
region has positive value. For Aef (r), the blue region has negative value (i.e. stabilized 
region) and the red region has positive value (i.e. destabilized region). Each panel shows 
the B^O distance of (a) 4.8 A, (b) 4.2 A, (c) 3.4 A, (d) 2.7 A, (e) 2.2 A and (f) 1.693 A 
(RC). 
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Figure 7: Same as Fig. [6] for the process from CI3B — OH 2 (RC) to the HC1 detachment. 
We partition the system into BC^OH) and HC1 for calculating An(f) and Aef(f). Each 
panel (a)-(g) corresponds respectively to the step 1-7 indicated in Fig. [ij The panel (c) 
describes TS. 
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Reaction steps 

Figure 8: The pathway of reaction of Fe(OH) 3 with BC1 3 toward Fe(Cl)(OBCl 2 ) 2 (OHBCl 2 ) 
(Sio) andFe(Cl) 2 (OBCl 2 )(OHBCl 2 ) (S 14 ). 
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Figure 9: The reaction pathways (a) to Fe(OH) 2 (OBCl 2 ) • HC1 (S 3 ) and to 
FeCl(OH) 2 (OHBCl 2 ) (S 4 ), (b) to Fe(OH)(Cl)(OBCl 2 )(OHBCl 2 ) (S 7 ), (c) to 
Fe(Cl)(OBCl 2 ) 2 (OHBCl 2 ) (Si ), and (d) to Fe(Cl) 2 (OBCl 2 )(OHBCl 2 ) (S 14 ). 
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Figure 11: Same as Fig. [6] for the process S5 — > TS2b —> S6 (see also Fig. [9] (b)). We 
partition the system into Fe(Cl)(OH)(BCl 2 0)(BCl 2 OH) and HC1 for calculating An(r) 
and Aef(f). TS2b is denoted by the panel (c). 
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Figure 12: Same as Fig. [6] for the process S2 — ► TSi —> S3 (see also Fig. [9] (a)). We 
partition the system into Fe(OH) 2 (BCl20) and HC1 for calculating An(f) and Aef(r). 
TSi is denoted by the panel (c). 
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